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We calculate the self-consistent in-medium T matrix for symmetric nu- 
clear matter using realistic interactions with many partial waves. We find 
for the interactions used (CDBonn and Nijmegen) very similar results for 
on-shell quantities. The effective mass and the renormalization factor Zp 
at the Fermi momentum are given for a range of densities. 

PACS numbers: 21.65.+f 

The presence of a strong interaction between nucleons at short dis- 
tances requires the use of resummation methods for many-body calcula- 
tions in nuclear matter. Besides the Brueckner-Hartree-Fock (BHF) sum- 
mation of ladder diagrams |2| and advanced variational methods ||, |j , the 
self-consistent in-medium T-matrix approach became feasible in last years 
H i> 0) H q- The numerical method presented in || allows for efficient 
calculations using off-shell nucleon propagators with realistic interactions. 



Below we present results for the CDBonn |L(| and Nijmegen potentials [11 
incorporating all partial waves with total angular momentum J < 9. The 
in medium T matrix fl 



T = V + VGGT (1) 

is calculated with the Green's function G = (u> — p 2 /2m — X) 1 self-consistently 
dressed by the self-energy in the T-matrix approximation 

iE = Tr[T A G] . (2) 

The difference with respect to the BHF approach lies in the dressing of the 
Green's functions by the full spectral functions, which means off-shell prop- 
agation. Also, the double Green's function in the T-matrix equation repre- 
sents the propagation of two particles or two-holes unlike in the Kernel of 
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the Bethe-Goldstone equation in the BHF scheme where the Pauli-blocking 
factor forces the propagation of two particles always. The above equations 
are solved by iteration for several nuclear densities p in the range 0.2 — 2.4 
normal nuclear densities po = .16fm -3 . The details of the equations and 
numerics can be found in Refs. |7|, |8|. 

We have improved the numerical algorithm to allow for many partial- 
waves in the T matrix. We use a separable parameterization of the inter- 
action, choosing 8 most important eigen-vectors of the interaction in the 
momentum representation for each uncoupled partial wave and 24 eigen- 
vectors for the coupled partial waves. This parameterization is essentially 
equivalent to the full parameterization in momentum. In previous T-matrix 
calculations a low-rank separable parameterization of the Paris potential 
was used. In this paper we compare the results for two realistic interactions 
without further simplifying approximations and taking many partial waves. 
With these calculations the binding energy and single-particle properties 
can be found without uncertainties due to technical simplifications. 

The binding energy per particle in the T-matrix approach can be calcu- 
lated from the Koltun's sum rule 

where A is the nontrivial spectral function obtained for the dressed prop- 
agators. In fact any expression for the energy should give the same result, 
since the self-consistent T-matrix approximation is thermodynamically con- 
sistent [13]. The existence a generating function $ for the self-energy guar- 



antees the fulfillment of thermodynamical relations between single-particle 
and global properties of the system. The fact that the T-matrix scheme is 
a consistent (conserving) approximation has been checked explicitely in 
using a model interaction. 

In the upper panels of Figs. [I] and § is shown the binding energy in the 
T-matrix approximation as a function of the density in symmetric nuclear 
matter compared to the corresponding BHF results. As noted previously 
the T matrix gives smaller binding energies and smaller saturation densities 
than the BHF calculation. At low densities the BHF and the T-matrix 
results converge as expected. The interactions studied in this work give 
similar results, e.g at normal nuclear density E/N = —14.3 and — 14.1MeV 
for the CDBonn and Nijmegen interaction respectively. The corresponding 
BHF binding energies at that density are 3.3 — 3.5MeV lower. In the same 
figures the Fermi energy Ep is shown. The Fermi energy in the T-matrix 
calculation is consistent with the binding energy; the Hugenholz-Van Hove 



relation [14| 



Ef = E/N (at saturation density ) (4) 
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Fig. 1. Upper panel : The binding energy and the Fermi energy as a function of 
density for the T-matrix approach and the binding energy in the BHF calculation, 
all for the CDBonn potential. Lower panel : The pressure as a function of density 
obtained from two different expressions Eqs. (||) and @. The solid and the dashed- 
dotted lines representing the two results for the pressure in the T-matrix calculation 
lie almost on top of each other. 



is automatically fulfilled. 

The pressure in the system can be obtained from several equivalent 
expressions. We consider 

(5, 

and 

P = p (Ep — E/N) . (6) 
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Fig. 2. Same as Fig. |T] but for the Nijmegen potential 



From the second form follows the Hugenholz-Van Hove relation at the sat- 
uration point (P = 0). The expressions (|5|) and ([]) are equivalent in the 
self-consistent T-matrix approximation. In the BHF approach the pressure 
can be calculated as the derivative of the binding energy (||). 

One-body properties are determined by the self-energy. In Fig. |3] is 
shown the imaginary part of the self-energy off-shell — Im£(p, uj). The self- 
energy is very similar for the CDBonn and Nijmegen potentials and close to 
the result for the Paris potential [|| for energies close to the Fermi energy 
(|w| < 200MeV) and for momenta up to 500MeV. At higher energies differ- 
ences start to appear, because of different off-shell behavior of the T-matrix 
for different interactions. The Paris interaction we use is a separable param- 
eterization with a small number of partial waves, this can lead to additional 
differences as compared to results from the other two potentials. 
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Fig. 3. The imaginary part of the retarded self-energy — Im£(p, lu) as a function of 
the energy oj for p — (upper panel) and p = 340MeV (lower panel) for different 
nucleoli- nucleon potentials. 



As can be seen from Fig. |I| the self-energy on-shell, i. e. for uj = u> p = 
p 2 /2m + Re£(p, uj v ), is similar for different interactions. For momenta p < 
700MeV the width of the quasiparticle excitation is similar. It is always zero 
at the Fermi momentum and increases quadratically when going away from 
it. The quadratic increase of the single-particle width is given by the cross 
section and the density of states at the Fermi surface. The cross sections 
are similar for different parameterizations of the nuclear potentials and the 
density of states at the Fermi surface is determined by the effective mass 
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Fig. 4. The imaginary part of the retarded self-energy — ImS(p, lo p ) at the quasi- 
particle pole as a function of the momentum. 



and the renormalization factor 

Z p = (l- dRe£(p, u)/du) | w=WjJ . (7) 

As shown below, those quantities obtained for different nucleon-nucleon 
interactions are similar. 

The real part of the self-energy determines the single-particle potential, 
i.e. the real part of the optical potential. The fulfillment of the Hugenholz- 
Van Hove relation guarantees the correct normalization of the single-particle 
potential. The real part of the self-energy is given as the sum of the Hartree- 
Fock and dispersive contributions 

R fiE (p )W ) = E OT (p)+/^^fc^. (8) 

J 7T UJ — U 

The Hartree-Fock part is different for the Nijmegen and CDBonn potentials, 
but the total single-particle energy on-shell is very similar (Fig. ||). The 
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Fig. 5. The single-particle potential Re£(p, u> p ) at the quasiparticle pole as a func- 
tion of the momentum. The upper curves represent the Hartree-Fock contribution 
for the Nijmegen and CDBonn interactions. 



difference in the Hartree-Fock part is compensated by the dispersive part, 
which includes integration of the imaginary part of the self-energy far off- 
shell. The differences in Im£(p, u>) for u> > 200MeV (Fig. ^) give the 
necessary shifts leading to similar total optical potentials (||). Although the 
Nijmegen and CDBonn interactions show some differences in the off-shell 
behavior of the self-energies, the resulting properties of the quasiparticle 
pole are very similar in the range of densities p < 2p$. 

The properties of the quasiparticles at the Fermi surface are defined by 
the renormalization factor (|7|) and the effective mass 

* P d P fn\ 
m = - — . (9) 

dujp 

The effective mass and the renormalization factor are important to define 



the effective interaction between quasiparticles ITR 16] and the superfluid 



gap [17, pTq] . We present the calculation of these quantities for different 
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Fig. 6. The effective mass m* /m at the Fermi surface for different nuclear interac- 
tions. 



interactions in a range of densities. We find that for p < 2po the CDBonn 
and the Nijmegen interactions give very similar results (Figs. |(| |7|). The 
difference of the calculation using the separable Paris potential can be at- 
tributed to the simplicity of that parameterization. The effective mass in 
our calculation comes out very close to the free nucleon mass. It is useful 
for further applications to parameterize its density dependence 

* / \ 2 

— = 1.43 -.623-^- + .146 f (10) 
m po \poJ 

and similarly the for renormalization Z p at the Fermi momentum (Fig. 

2 



Zp = .57 + .2— — .05 f — I . ( I I) 



At normal nuclear density Zp ~ 0.72 which means a reduction of the effec- 
tive interaction between quasiparticles by a factor Zp ~ 0.5. 
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Fig. 7. The renormalization factor of the quasiparticle pole Z v at the Fermi surface 
for different nuclear interactions. 



We extend previous calculations of the self-consistent T matrix to soft 
core CDBonn and Nijmegen nucleon-nucleon potentials. We include in the 
calculation partial waves up to the total angular momentum J = 8. We 
present for the first time results on the binding energy for such realistic and 
detailed parameterizations of the two-body interactions. The binding energy 
leads to a harder equation of state, a smaller binding energy and a smaller 
saturation density than the BHF approximation. The pressure we find is 
thermodynamically consistent as expected for a conserving approximation. 
The single-particle energy and the width at the quasiparticle pole come out 
similarly for different interactions used. It is remarkable that differences 
in the Hartree-Fock energies and differences in the off-shell behavior of the 
imaginary part of the self-energy cancel out in the single-particle potential. 
We present results on the properties of the quasiparticle pole at the Fermi 
surface. The effective mass is close to the free one and the renormalization 
factor of the quasiparticle pole is Zp ~ 0.7 around the normal nuclear 
density. 
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